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Abstract 


By  means  of  a  sampling  theorem,  the  passage  is  made  from 
lattice  dynamics  to  nonlocal  elasticity.  Foundations  of  nonlinear  (and 
linear)  nonlocal  elasticity  is  established  on  the  basis  of  the  nonlocal 
actions  in  continuum  theory.  The  solutions  are  presented  for  the  pro¬ 
blems  of  screw-dislocation,  anti-plane  line  crack  and  point  defect  in  an 
elastic  solid.  Agreements  with  the  atomic  results  and  experiments  are 
remarkable.  -  j 


Introduction 


All  materials  are  made  up  of  subbodies  which  constitute  their  sub¬ 
structure  or  microstructure.  The  size  of  subbodies  may  vary  from  atomic 
dimensions  to  the  macroscopic  scale  such  as  grain  size.  Depending  upon 
the  nature  and  accuracy  of  the  physical  phenomena  to  be  modelled,  the 
average  distance  of  subbodies  play  a  central  role.  This  distance  may 
vary  from  the  order  of  the  lattice  parameter  a  (10“8  cm  in  perfect  cry¬ 
stals),  to  a  few  millimeters  as  in  granular  solids.  The  boundary  and 
initial  conditions  brings  into  play  another  characteristic  length  A 
(e.g.  thickness  of  a  plate,  wave  length).  The  domain  of  applicability 
of  a  mathematical  model  depends  on  the  ratio  e  = a/A.  For  classical 
continuum  theory  e  =  0  and  for  the  lattice  dynamics  e  =  1  .  Since  for 
real  amorphous  materials,  lattice  dynamic  calculations  are  not  possible 
or  extremely  costly,  we  may  ask  whether  continuum  theories  can  be  con¬ 
structed  to  cover  the  range  0<  e<  1. 


Indeed,  such  models  entitled  micropolar,  micromorphic  and  nonlocal 
theories,  have  been  proposed  during  the  past  two  decades.  Literature  is 
fairly  extensive  on  several  sections  of  these  theories  (see  [1,2]).  The 
solutions  of  various  critical  problems  have  verified  our  hopes  and  ex¬ 
pectations  in  that  by  means  of  nonlocal  models,  accurate  predictions  are 


possible  of  the  physical  phenomena  in  submicroscopic  scales.  As  a 
result,  various  long  standing  problems,  such  as  crack  tip  singularity, 
dislocation  core,  point  defects  have  been  treated  accurately  [3-6], 
Nonlocal  plasticity  [7] ,  nonlocal  fluid  mechanics  [8] ,  nonlocal  electro¬ 
magnetism  [9,10]  are  yet  to  be  explored.  Nevertheless,  a  few  problems 
treated  in  these  areas  have  already  shown  their  power  and  potential, 
indicating  that  we  are  on  tlie  right  path. 

In  this  paper,  first  I  present  an  account  on  linear  nonlocal  elas¬ 
ticity  based  on  lattice  dynamical  models  (Section  2) .  This  shows  us 
the  firm  physical  foundations  of  the  nonlocal  elasticity.  Next,  I  pre¬ 
sent  a  general  theory  of  constitutive  equations  (Sections  3  and  4) , 
which  can  be  used  to  develop  constitutive  equations  of  memory -dependent 
materials  (e.g.  polymeric  substances)  undergoing  finite  deformations. 
Section  5  contains  an  account  of  some  approximate  models  which  can 
account  for  shurt  nonlocality.  These  models  are  conducive  to  simpler 
mathematical  treatments.  The  last  three  sections  (Sections  6-8),  con¬ 
tain  solutions  of  some  critical  problems  which  fall  outside  of  the 
domain  of  classical  elasticity,  since  they  give  rise  to  stress  singu¬ 
larities.  The  solution  of  the  screw-dislocation  problem  (Section  6) 
leads  to  a  finite  maximum  stress  near  the  dislocation.  The  cohesive 
strength  of  perfect  crystals  based  on  this  solution  is  in  remarkably 
good  agreement  with  those  based  on  the  atomic  theories.  The  solution 
of  the  anti-planf  crack  problem  presented  in  Section  7  contains  no 
stress  singularity  so  that  the  fracture  criterium  can  be  based  on  phy¬ 
sically  realistic  maximum  stress  hypothesis.  As  shown  in  our  other 
works  (cf.  [3,5,11]),  it  gives  excellent  agreement  with  experimental 
results.  The  last  section  contains  the  solution  of  the  point  defect 
problem. 

2.  Lattice  Dynamical  Foundation  of  Linear  Elasticity 

According  to  lattice  dynamics,  a  perfect  crystal  consists  of  discrete 
atomic  mass  points  attached  to  each  other  by  springs  (Fig.  la).  Let 
u(n)  denote  the  displacement  vector  of  an  atom  located  at  a  discrete 
point  marked  by  n  from  the  equilibrium  position  x°(n) .  If  e^  are 
a  set  of  base  vectors  and  ek  are  their  reciprocals,  i.e. 

Ik  o  _  xk 

2  *  ~  6  l 


We  may  write  for  the  undeformed  and  deformed  positions  of  an  atom 
initially  located  at  x°(n)  (Fig.  lb); 


(2.1) 


2k£k  ’ 


x(n)  =  x  (n)  +  u(n) 


For  a  net  consisting  of  atoms  which  are  located  at  the  "knots"  of 
a  parallel  piped  with  equal  edge  lengths  (which  are  parallel  to  e,) 
nK  may  be  taken  to  be  integers.  Such  a  lattice  is  called  a  simple 


The  potential  energy  (internal  energy)  of  a  lattice  is  a  function 
only  of  the  distance  between  atoms,  e.g. 


■  FSM..TZM  iGLJUSVJUiai  '&a*JVK*Vm.RZM  _r 


(2.2)  v  =  l  U[jx(n)  -  x(n»)|] 

'  '  |  *w  *1/  «v 

n  n* 

«V 

A  Mechanical  analogy  is  a  spring  where  the  restoring  force  is  a  func 
tion  of  the  change  of  length  of  the  spring.  For  the  harmonic  approxi 
mat ion,  it  is  simple  to  show  that  only  quadratic  terms  in  U  ,  in  the 
polynomial  expansions  of  U  about  u  =  0  survives,  and  we  have 


uk(n)  $u(n,n’)  u£(n’) 


(2.3) 


V 


where  the  constant  potential  energy  corresponding  to  the  equilibrium 
configuration  is  discarded  and 


(2.4)  \£(n,n:)  =  U>kJl(|x(n)  -  x0(n*)  |) 


The  Lagrangian  of  the  system  has  the  form 

(2.5)  L  =  i  l  m(n)u(n)«u(n)  -  i  I  u(n>  $(n,n’)«  u(n’) 
zn .  ^  nn'  ~  ~  ~  ~  ~ 

-  I  q(n)«u(n) 

ni  >v 

n 


where  m(n)  is  the  mass  of  ncn  atom  and  o(n)  is  the  external  force 
applied  to  it. 

Lagrange’s  equation  now  reads 

(2.6)  m(n)  u  +  l  4>(n,n’)*u(n’)  =  q(n) 

<v  »w  «v 

n’  ~ 

This  is  the  basic  equation  underlying  the  non-relativistic  lattice 
dynamics. 

To  pass  to  the  so-called  quasicontinuum  representation,  we  need  to 
convert  the  summations  in  (2.5)  and  (2.6)  to  integrations.  This  is 
done  by  means  of  a  sampling  function  based  on  a  sampling  theorem  in¬ 
troduced  by_Shannon  [12 J .  According  to  this  theorem,  if  the  Fourier 
transform  f(k)  of  a  function  f(x)  vanishes  outside  a  region 
|k|  >  k  ,  then  this  function  is  determined  uniquely  from  its  discrete 
values  f(n)=f(x)  at  x  =  n  .  In  fact,  if 

(2.7)  f(k)  =  v0  l  f(n)  e^’lS 


is  the  Fourier  transform  of  f(x)  ,  then  the  inverse  transform  is  given 
by 

(2.8)  f(x)  =  vQ  l  f(n)  6B(x-n)  , 

2 

where 

(2.9)  <SR (x)  =  -i--  e%2  dk 

a  ~  (2TT)-5 

B 

and  v  is  the  volume  of  an  elementary  cell  constructed  on  the  vector 
e^  .  For  a  simple  lattice  with  lattice  parameter  a  ,  |e^|  =  a. 

v  =  a3  and  the  domain  of  k  is  0<  |k|<  1.  It  is  clear  that  6fi(0)  = 

v0  *  and  6R(n)  =0  for  all  other  n^O. 

D  <v  +* 
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We  note  the  important  fact  that  f(x)  given  by  (2.8)  coincides 
exactly  with  f(n)  at  discrete  points  x  -n,  however  being  completely 
arbitrary  elsewhere.  Thus,  we  expect  that  tihe  predictions  of  the  non¬ 
local  theory  will  agree  with  those  of  lattice  dynamics  at  the  atomic 
points,  differing  possibly  in  between. 

Using  the  foregoing  process,  we  obtain  the  nonlocal  representation 
of  (2. 3) and  (2.6) 

(2.10)  V  =  i  u(x)  •  <Kx,x,)*u(x')  dv(x)  dv(x') 

4  A«  M  ^  A#  <W  Ai  ^ 
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(2.11)  pii(x)  +  ‘KxjX’^utx’)  dv(x’)  =  q(x) 


where  p(x)  is  the  mass  density. 

The  invariance  of  $  under  ri^.id  translations  and  rotations 
further  shows  tnat  (cf„  Kunin  [  6,  p.  14],  the  Fourier  transform  of 
4  can  be  written  as: 


(2.12) 


<Tj  (k,k»)  =  (k,k’)  . 


Consequently, 

(2.13) 


V  = 


c  Cx,x’)  eu(x)  emj(x’)  dv  dv»  , 


(2.14) 


t  v  +  p(r  -  u  )  =  o 

»x 


where  we  wrote  dv' =  dv(x')  and 


(2.15) 


(2.16) 


£  _  -£ 
q  -  p£  » 


emn  1  ^um,n  +  un,m^ 


Here  represents  functional  partial  derivative.  If  the  ma¬ 

terial  properties  vary  slowly  over  a  distance  d  ,  then 

(2.17)  ck!“  C2,x')  =  ck!“(x)  SgCx-x1) 


5 


it*  icnjuT  a  t^-.* 


where  cl^Hx)  are  the  elastic  moduli  of  zero1"  -order,  with  respect 
to  d  .  In  the  limit  as  the  internal  characteristic  length  (e.g. 

a)  goes  to  zero,  SR  becomes  Dirac  delta  measure  and  (2.15)  gives 
classical  Hooke’s  law.  Now  we  are  on  familiar  grounds. 

The  function  6g(x)  given  by  (209  )  have  been  determined  for 
various  crystals.  For  cubic  lattices,  6g  has  the  simple  form 


(2.18)  6g (x-x ’ ) 


i=l  xrx’ 


sin 


a 


Another  simple  case  is  that  of  Debye  continuum,  which  corresponds  to 
the  case  of  a  spherical  Brillouin  zone,  with  radius  <  =  m/a  .  For  this 
case,  we  have  (cf.  Kunin  [6,  p.  37]) 


(2.19)  6B(x)  =  -^y~2  -  cos  <r)  , 

2  m  r 


r  ~  v6c*x 


for  the  three-dimensional  case  and 

(2.2°)  5B(x)  =  ^  Jx(Kr) 


for  the  two-dimensional  case. 


For  engineering  materials,  it  is  useful  to  select  other  functions 
which  give  better  match  with  the  dispersion  curves  in  the  neighborhood 
of  k  =  0  .  We  make  the  following  important  observations : 


i) 

ii) 

iii) 

iv) 

v) 


In  nonlocal  elasticity,  the  total  internal  energy  is  a 
functional  of  the  strain  (eq.  2.13) . 


The  internal  energy  density  coincides  with  those  of  the 
lattice  dynamics  at  the  discrete  points  occupied  by  the  atoms. 


The  displacement  fields  at  the  atomic  points  are  the  same 
for  both  the  nonlocal  elasticity  and  the  lattice  dynamics. 


Nonlocal  elastic  moduli  can  be  determined  from  the  inter¬ 
atomic  potential. 


In  the  limit,  as  the  internal  characteristic  length  (e.g. 
lattice  parameter  a )  goes  to  zero,  nonlocal  elasticity 
reduces  to  classical  (local)  elasticity. 


The  development  given  above  is  for  the  so-called  haimonic  approximation 
(linear  theory)  and  non-dissipative  systems. 


The  interatomic  potential  is  a  nonlinear  function  of  the  distance. 
The  harmonic  approximation  is  valid  only  in  the  neighborhood  of  its 
minimum.  When  dislocations,  holes  and  impurities  exist  (as  in  the  case 


i,t.  .r  «?r . ;r.  ^  -r.  ,T.y,  •_  v.v  v  v -,-y 


of  real  materials) ,  the  dissipative  effects  come  into  play.  In  such 
situations,  it  is  very  difficult  (if  not  impossible)  to  construct  a 
nonlocal  theory  starting  from  lattice  dynamics.  Moreover,  continuum 
approach  make  sense  on  its  own  grounds  as  an  independent  discipline. 
In  Section  4,  we  present  one  such  development. 


Basic  Axioms  of  Constitutive  Theoi 


Some  years  ago  (cf.  Eringen  [13]),  I  organized  the  basic  elements 
of  the  constitutive  theory  of  continuum  mechanics  into  eight  axioms: 


(i)  Axiom  of  Causality 

(ii)  Axiom  of  Determinism 

(iii)  Axiom  of  Equipresence 

(iv)  Axiom  of  Objectivity 

(v)  Axiom  of  Material  Invariance 

(vi)  Axiom  of  Neighborhood 

(vii)  Axiom  of  Memory 
(viii)  Axiom  of  Admissibility 


Of  these  (vi)  and  (vii)  represents  approximations  that  can  be  made  be¬ 
cause  of  the  decay  of  interatomic  interactions  with  distance  and  past 
times.  Axiom  (i)  enunciate  that  the  motion  and  temperatures  of  all 
material  points  are  the  cause  of  all  physical  phenomena.  In  fact,  even 
the  temperature  can  be  discarded  if  the  intrinsic  motions  of  sub- 
bodies  are  taken  into  account.  For  example,  according  to  the  kinetic 
theory,  the  average  kinetic  energy  of  intrinsic  motions  is  the  absolute 
temperature. 


Axiom  (ii)  and  (iii)  state  that  all  other  state  variables  (e.g. 
internal  energy) ,  at  a  point  of  the  body,  are  functionals  of  the  motions 
and  temperatures  of  all  points  of  the  body.  Here  is  the  nonlocal idy ! 

Axiom  (iv)  and  (v)  expresses  the  invariance  requirements  for  the 
response  functionals,  the  first  being  the  invariance  under  translations 
and  rotations  of  the  spatial  frame  and  the  second,  under  the  material 
symmetry  group. 

Finally,  Axiom  (viii)  imposes  the  compatibility  of  these  func¬ 
tionals  with  the  fundamental  laws  of  motion  and  thermodynamics. 

In  Section  4,  we  proceed  to  show  how  these  ideas  can  be  used  to 
develop  a  nonlocal  continuum  theory. 


of  Nonlocal  Elasticity 


According  to  the  axiom  of  causality,  all  state  variables  Z  (Free 
energy  \p  ,  entropy  q  ,  stress  ly,  ,  heat  %)  are  functionals  of  the 
motions  and  temperatures  of  all  points  X’  of  the  body,  at  all  past 
times  - 00  <t’  <  t  .  For  simplicity,  here  we  consider  only  isothermal 
elastic  solids  with  no  memory.  Consequently,  we  have 


(4.1)  *(X,t)  =  ¥[x(X',t),X*]  , 


X’  el/ 


V¥W. 


where  x’  =  x(X',t)  is  the  motion  of  a  point  X’  in  the  body  V,  and 
¥  is  a  functional  of  the  motion  of  all  points  X*el/  .  This  is  also 
postulated  for  all  other  members  of  Z  .  Applications  of  the  invariance 
under  rigid  body  translations  and  rotations  shows  that  ^  is  a  func¬ 
tional  of  the  distance  between  the  reference  point  and  all  other  points 
in  the  body 


(4.2) 


<KX,t)  -  ¥[|x(X*,t)  -  x(X,t)  | ,  X] 


But  we  have,  in  rectangular  coordinates  (see  [13,  p.  45]), 


(4.3)  xk(X’,t)  -xk(X,t)  =  j  R^C^dXK 

X 


where  Rj,  and  are,  respectively,  the  rotation  and  defomiation 

tensors.  Consequently,  tp  can  be  considered  a  functional  of  g  an.i 
C  .  The  dependence  on  R  is  admissible  only  through  the  difference 
R(X' ,t)  -  R(X,t)  on  account  of  the  objectivity.  For  non-polar  ma¬ 
terials  "(Axiom  (v)) ,  this  dependence  is  eliminated  completely  so  tliat 
we  arrive  at 

(4.4)  ip(X,t)  =  *[C(X»,t),X*  ;C(X,t),X]  ,  X7X 


For  non-heat  conducting  solids,  the  second  law  of  thermodynamics 
states  that 


(4  c  5) 


[-  Po(i  +  ri0)  +  tklclk]  dv  >  o 


where  pQ  ,  n  ,  0  >  0  and  T^  are,  respectively,  the  mass  density  of 

the  body  at  the  natural  state,  the  entropy  density,  the  absolute  tem¬ 
perature  ,  and  the  Piola-Kirckhoff  stress  tensor  of  which  T^.  is  rela¬ 
ted  to  Cauchy  stress  tensor  tk£  by  ^ 

(4o6;)  tldt  =  p“  TKL  xk,Kx£,L 


Note  that,  as  against  classical  (local)  continuum  theory,  we  posit  that 
the.  total  entropy  of  the  body  must  be  non-negative. 

The  axiom  admissability  (viii)  requires  that  (4.4)  must  not  violate 

(4.5).  In  order  to  calculate  <J>  we  need  some  smoothness  requirement 
for  the  functional  T.  To  this  end,  we  can  assume  that  C(X*,t)  belongs  to 
a  Banach  or  Hilbert  space.  However,  if  we  confine  our  appe¬ 

tite  to  a  more  limited  class  of  functionals,  namely  the  linear  func¬ 
tionals,  in  the  sense  of  Friedman  and  Katz  [14],  then  we  can  write 


(4.7)  pto(X,t)  =  f  F[C(X*  ,t),X*  ;C(X,t),X]  dV(X’) 

(J  **  Mm  w  <s* 


The  total  internal  energy  is  given  by 


(4=8) 


pKX,t)  dV 


F [C *  ,X'  ;C,X]  dV'  dV 

*V  N  M 


V  V  V 

Since  the  change  of  the  order  of  integrations  does  not  affect  V ,  we 
assume  that 

(4.9)  F[C’,X’  ;  C,X]  =  F  ,  F  =  F[C,X;C\X'] 

*v  iv  »v  ^  <y 

where  and  henceforth,  a  superposed  asterisk  indicates  the  interchange  of 
X  and  X'  . 

We~assume  that  F  is  continuously  differentiable  with  respect  to 
C*  and  Co  Using  (4.7),  we  calculate 


If  this  is  to  be  valid  for  all  values  of  C  throughout  the  body, 
it  follows  that 

('4°14')  Tkl  =  +  ^4^  ^  dV’ 

Thus,  we  have  proved 

THEOREM:  The  stress  constitutive  equations  given  by  Equation 
(4. 14)  is  compatible  with  the  second  law  of  thermo¬ 
dynamics. 

Note  that  (4014)  is  exact,  and  it  is  valid  for  arbitrarily  large  elastic 


imm  ’if*  SflCtTVi^ELHI^Vi^ri  .  xjwwv: 


■*/mi  wv\  wt^u* 


strains .  By  means  of  polynomial  expansions  of  F  in  C'  ,  we  obtain 
various  nonlinear  theories  and  the  linear  theory. 


Linear  Theory 

It  is  practical  to  introduce  the  linear  strain  measures 


(4,15) 


CKL  *  SKL  +  2  EKL  ’ 


EKL  '  7  'UK,L  *  \,V? 


where  U  is  the  displacement  vector.  We  write 

<4‘16>  F  -  co  + Wkl% 

where 


(4.17) 


C0  =  C0  ScLMN 


~  CTT/.At  =  c 


LKMN  KLNM 


Sm 


Substituting  (4.16)  into  (4.14),  we  obtain 


(4,18) 


KL 


ww)  w>«  «■ 


For  the  linear  theory,  (4,6)  gives  hi  ~  TKL  6kX6  ,  consequently, 


(4.19) 


hi 


ero(i'>t)  dvfx') 


where 

(4.20) 


:k£mn  ”  CKLMN  6Kk6L&VNn  » 


1  ,  ,  s 

1-U  t,  +  J  . 


“mri  7  ^  m,n  n,nr 


Hence  we  obtained  the  constitutive  equations  of  the  linear  theory  of 
nonlocal  elasticity  again,  without  any  resort  to  the  atomic  structure. 


For  homogeneous  medium  c  =  c(x'-x)  and  for  the  homogeneous  and 
isotropic  medium,  c  will  be~an~isotropic  function  of  k  =  x*-x  ,  i.e. 

**»  •#  M 


(4.21) 


"kEinn  ~  Xo  dk;d'mn  +  ^o^km^En  +  dkndEm-* 

‘Vw  {k«  *  ViW 


+  VVm  V  +  "k^M  *  Vm5kn+ 


+  Vn  V 


+  x3  KkKE  Vn  > 
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where  X  .y„  and  X  are  functions  of  I k I 
o  ’  o  a  1 «. 1 

(4.22)  (X0,P0>X0)  =  U0(|£|),  V>0 C I k I ) »  X0(||C|)  , 


a=l,2,3 


The  material  moduli  Xa  are  the  contributions  of  the  interatomic  orien¬ 
tations  to  the  isotropy.  They  are  not  present  in  local  elasticity.  Be¬ 
cause  of  the  appearance  of  it,  we  expect  that  these  terms  are  gener¬ 
ally  much  smaller  than  XQ  aria  yQ  .  We  also  expect  that  XQ  and  yQ 
attenuate  with  distance  in  the  same  way.  If  these  two  assumptions  are 
accepted,  then  (4.19)  may  be  put  into  the  simple  form 


(4.23) 


where 


’ 

a(|x'-x|)  ok£(x')  dv(x’) 

1/ 


is  the  classical  (local)  Hookean  stress,  i.e., 


(4.24) 


aki  =  x  err<ska  +  2veU 


and  X  and  y  are  the  classical  Lame  constants. 

The  kernel  function  a(|x'-x|)  appearing  in  (4.23)  is  of  the  di¬ 
mension  length" 3.  Consequently," it  must  depend  on  the  internal  charac 
teristic  length.  In  fact,  according  to  Eq.  (2.17), 

(4.25)  cx(|x’-x|)  =  6b(x'-x) 


5.  Approximate  Models 

(i)  If  the  fields  are  slowly  varying  over  a  characteristic  dis¬ 
tance  d  ,  we  can  replace  the  integral  operators  with  differential 
operators.  This  simply  means  that  we  expand  c(k)  or  a(k)  into  power 
series  of  k-  ,  e.g. 


(5.1) 


ck£mnT.,  ...t„ 
1  r 


(i  kj)  (ik2) . . . (ikp) 


With  ck£mnT  being  real  constants,  (4.19)  gives  couple  stress 

1’  *  *  r 

theory  involving  r  derivatives  of  e^  . 


(ii)  A  more  useful  case  involves  matching  c  or  a  in  the  wave 
number  space  with  the  dispersion  curves,  based  on~atomic  models.  For 
example,  according  to  the  Bom-Karman  lattice  model, 


(5.2) 

where 

(5.3) 


w?(k)/w^.  =  (2</irk)^  sin^(iTk/2K) 

4>Oj  *  Af  >  -  O+Wp  ,  z\  =  p/p 


are,  respectively,  the  classical  circular  frequencies  of  the  irrotation 
al  and  equivoluminal  waves.  Here  k  is  the  upper  limit  of  k  (the 
boundary  of  the  Brillouin  zone),  e.g.,  for  one-dimensional  lattice 
k  =  ir/a  ,  where  a  is  the  lattice  parametei . 

If  we  use  (4,23)  with  (4.24)  in  (2.14),  with  f  =  0  for  the  har¬ 
monic  waves,  we  obtain 

(5.4)  w?(k)/e£.  -  a(k2) 

where  a  is  the  Fourier  transform  of  a  .  A  polynomial  form  of  1/a 
is  most  convenient  to  approximate  (5.2): 

(5.5)  1/ci  =  1  +  3j(k/«:)2  +  (k/<)4 


At  k  =  0  ,  this  gives  classical  value  =  o>0 j  and  at  k  -  <  ,  it 
satisfies  the  condition  dw(k)/dk  =  0  if  6.=1.  Hence,  we  have 


(5.6) 


(5.7) 


S  -  [1  +  SjCk/K)2  +  (k/K)4]"1 

«/uoj  •  [1  *  Bjfk/ic)2  *  (k/K)4]"1/2 


The  parameter  $,  may  be  connected  with  the  ratio  of  the  boundary 
frequency  w. (k)  to  J  the  Debye  frequency  w4  ,  i.e. 

4  1 


(5.8) 


WjGO/Wj  =  (2  +  Bj)”^2 


from  which  it  follows  that  -2  <  g.  <®  .  pig.  2  displays  the  dispersion 
curve  based  on  (5.7),  along  with  J  that  based  on  the  Bom-KSrmSn  model 

V.  J  ^  J  •  j 


♦-BORN-KARMAN 


0  0.2  0.4  0.6  0.8  1.0 

k/K 


Figure  2:  Dispersion  Curves 


AfUKjnua  r. ucx  1. u  .x  u 


iwr -»  *r- 


It  is  more  practical  to  leave  3*  and  6j  free  for  a  better  fit¬ 
ting  at  a  lower  wave  numbers  for  macroscopic  engineering  problems. 

Therefore,  we  take 

(5.9)  a  =  (1  +  e2k2  +  yV)"1  . 

With  this  we  have 

(5.10)  (1  +  e2k2  +  YV)  tk£  =  ok£ 

The  inverse  Fourier  transform  of  this  gives 
(5.U)  (1  -  e2V2  +  yV)tk£  =  on 

In  terms  of  (4.23),  this  implies  that 

(5.12)  (1  -  e2V2  +  y^V*)  a(x'-x)  =  6(x*-x) 

*V  *W  A# 

where  S(x’-x)  is  the  Dirac  delta  function.  If  we  take  the  divergence 
of  (5.11)~an3  use  (2.14),  we  obtain 

(5.13)  ak£jk  +  (1  -  e2V2  +  y4V4)(p  f£  -  pii£)  =  0 
or,  on  using  (4.24), 

(5.14)  (X  +  2u)  +  U  ul(kk  ♦  (1  -  e2V2  +  yVXpf*-  pip  -  0 


These  equations  replace  the  Navier  equations  of  classical  elasticity. 
They  are  singularly  perturbed.  Note  that  for  the  static  problems  and 
vanishing  body  forces  (5.13) reduces  to  the  equations  of  equilibrium 
of  the  classical  elasticity  and  (5.14)  to  Navier* s  equations.  However, 
since  the  real  stress  is  not  Oy.  but  tk£,  in  order  to  determine  the 
stress  field,  we  must  invert  (5.10). 

6.  Screw  Dislocation 

A  screw  dislocation  is  a  constant  discontinuity  in  the  displace¬ 
ment  component  u^  at  a  point  in  the  plane  x, =  0  .  The  discontinuity 
vector  has  a  single  component  of  magnitude  b  nailed.  Burger’s  vector. 
In  this  case,  u^=U2  =  0  and  u?  =u,(x.,  ,x2)  .  Since  all  components  of 
the  stress  tensor  vanishes  except  t^,  and  t_2  ,  the  equations  of  equi¬ 
librium  are  satisfied  by  introducing’5 the  stress  function  4>(XpX2)  by 

(6.1)  *13  =  ^,2  ’  ^23  ^1 

Equations  (5.11),  in  the  simplest  case  y  =  0,  gives 

(6.2)  (1  -  e2V2H>2  =  o13  ,  (1  -  eV)*^  -  - 
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we  have  the  expressions 


°13  ‘  >,u3,l- 


ff23  ‘ 


For  o13  and  a^, 

(6.3) 

from  which  it  follows  that 

(6.4)  a 

Using  (6.2),  this  gives 

(6.5)  (1  -  e2v2)  V2<j)  =  0 


3,2 


13,2  a23,l 


An  appropriate  solution  of  (6.5)  vanishing  as  r=  (xf  +  x?)^  ->  oo  is 
given  by  11 

C6-6)  <P  =  CQ  +  C^n  r  +  C2  K0(r/c) 

where  Kq (z)  is  the  modified  Bessel’s  function,  Cq,  and  C2  are 
arbitrary  constants,  and  (r,6)  are  the  plane  polar  coordinates,  i. 


(6.7) 


x^  =  r  cos  0  ,  x2  =  r  sin  0 


The  stress  field  referred  to  polar  coordinates  is  given  by 
(6.8)  tzr  *  t31  cos  0  +  t32  sin  0, 

tzQ  =  -  t31  sin  0+  t32  cos  0 

It  is  single  to  show  that 


(6.9) 


_  1  3<J> 


f.  - 


34> 

9r 


zr  r  90  »  ''zd 

Employing  (6.2),  (6.3)  and  (6.6),  we  find  that 


(6.10) 


u3  =  ci 0 


The  jump  discontinuity  u3(2tt)  -  u3(0)  =  b  is  the  Burger’s  vector  b 
Consequently,  =  -b/2ir ,  and  we  nave 

(6.11)  <J>  =  Cq  -  (yb/2rr)  Jin  r  +  C2K^(r/e) 

At  the  eye  of  the  dislocation  r  =  0  ,  t  Q  will  vanish  if  C2  = 
-  yb/2ire  .  Hence,  the  stress  field  is  given  by 

<6-12)  tze  =  ft  [1  -  l^tr/e)]  ,  t2r  =  0 

This  result  was  also  found  in  another  way  in  our  previous  works  [5] . 


The  non-dimensional  stress  field 
(6.13)  T  =  (2-rre/b)  (tz0/u)  =  p**1  [1  -  pK^p)] 


is  plotted  as  a 
suit  t 

sical  zwvalue 
given  by 

(6.14) 


function  of  p  =  r/e  in  Fig.  3,  Unlike  classical  re- 
has  no  singularity  at  p  =  0  and  it  approaches  clas- 
1/p  for  large  p  .  It  has  a  maximum  at  about  p  =  l.l  , 


t,fi  /p  -  0.3993 
z0  max  2tte 


The  strain  energy  per  unit  length  L  in  the  z-direction  of  a  cylindri¬ 
cal  tube  of  inner  radius  r  and  outer  radius  R  ,  is  given  by 

R  0  9 

f  ..v2 

(6.15)  Z/L  =  2tt 


r  J 
Lo 


tze  ez6  T  dr  *  W-  [»0V«o)  +  K0(R/e)  -K0(ro/e)] 


Again,  contrary  to  the  classical  result,  this  has  no  singularity  at 
r0=  0  .  For  large  R  ,  the  strain  energy  depends  on  the  size  of  the 
solid,  as  expected. 

By  use  of  (6.14),  we  can  determine  the  cohesive  strength  t  of  a 
perfect  crystal.  To  this  end,  we  write 


(6.16) 


t  „  =  0.3993  =  t 

z0  max  2ne  c 


Based  on  phonon -dispersion,  e  -  0.39a,  where  a  is  the  lattice 

parameter.  Using  this  for  an  aluminum  crystal,  we  find  that 

(6.17)  tc/y  «  0.12  {M:  [111]  <110  >} 

This  is  very  close  to  the  theoretical  strength  ty/u  =  0.11,  based  on 
atomic  models.  For  distributed  dislocations  macroscopic  strength 
can  be  estimated.  For  such  considerations,  see  [15], 


7.  Mode  III  Crack 

A  line  crack  located  at  |x^j  <c,  x2  =  0  in  the  plane  x^O  is 

subject  to  a  constant  shear  stress  t23  =  -t0  alon2  UPP61  surface 
of  the  crack  and  tg  along  the  lower.  The  stress  field  is  that  of 
Mode  III,  i.e.  the  only  non-vanishing  stress  components  are  t15  and 
t23.  To  determine  the  stress  field,  we  must  obtain  appropriate  solu¬ 
tions  of  (6.5)  and  (6.3),  subject  to  the  boundary  conditions 


(7.1) 

t23  : 

u3  ! 

(7.2) 

{t23 

t0  » 


as 


ix1|  <  c  ,  x2  =  o 

|Xll  >  c  ,  x2  =  0 

(x*  +  x*)%  •>  00 
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Figure  3:  Non-Dimensional  Hoop  Stress  for  Screw 

Dislocation 


(7«3) 


*13  =0 


*l\  -’■c  , 


x2  =  ° 


Of  these,  (7.3)  expresses  the  vanishing  of  the  surface  traction  on  the 
tip  surfaces  as  these  surfaces  become  two  distinct  points  at  =  0  . 


An  appropriate  solution  of  (6.5)  vanishing  at  infinity  in  the 
upper  plane  X£  >_  0  ,  is  found  by  means  of  the  Fourier  transform  tech¬ 


nique 

(7.4) 


*  -  c|)S 


<x2 


-  Ax„ 


[?  A(0  e  z  +  B(0  e  *]  sinCSx^d? 


0 


where  A(g)  and  B(£)  are  to  be  determined  from  the  boundary  condition 
and 


(7.5) 


X  =  (E2  +  e'2)* 


Using  (6.2)  and  (6.3),  we  find  that 


00 


(7.6) 


The  stress  field  is  given  by 


u,  =  (3*  ZHZ)  e  cosfCx-,)  d£ 


*  **  Xx 

Wu  =  -  (f)*  [?2  A(0  e  2  +  AB(£)  e  2]  sin(£Xl)d£, 


(7.7) 


r^\h  (  rr2A,_-v  ^*2  .  enfO  ~  2-1 


t23/y  =  "  (£r  [rA(0  e  4  +  SB(£)  e  "]  cos^d? 


At  the  tips  of  the  crack,  t.^  vanishes  if 
(7.8)  B  =  A(0  Z2  A/X 


(7.8) 

and  we  obtain 


tjj/n  =  -  (§)*  }c2AK)Ce  ?*2  -  e  sinCCXj)  d(  , 


(7.9) 


2  J?  ?  ~Zi^2  p  “lXo 

t23/^J  =  -  C^r  J5AG)(e  -  f  e  2)  cos^)  d£ 


The  remaining  boundary  conditions  at  x2  =  0  read 


(7.10) 


where  we  set 


(7.11) 


w 

?  C(0  [1  +  k(eO]  J.*(xC)dC  =  t0/y  x% ,  0  <x  <  1 

* 

3 

00 

J  Cft)  J.%(x?)  d;  »  0  ,  X  >  1 


k(ec)  -  -  CU2  *  e2]"1/2 
x  =  x-j/c  ,  s  =  cZ 
C(?)  =  c"3  ?3/2  A(?/c) 


e  =  c/e 


and  used  Jt(z)  =  (2/ttz)^  cos  z  ,  to  replace  cos(x£)  with  its  expres¬ 
sion  in  terms  of  the  Bessel's  function  J_%(z). 

The  solution  of  the  dual  integral  equations  (7.10)  is  not  known. 
However,  it  is  possible  to  reduce  the  problem  to  the  solution  of  a 
Fredholm  equation  as  shown  in  our  previous  work  [11] .  An  approximate 
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solution  for  small  e  can  be  obtained  and  it  has  the  classical  form 

(7.12)  C0(C)  =  (if/2)%  (t0/y)  5*  ^(0 

This  solution  suffices  for  calculations  when  the  crack  size  is  larger 
than  100  atomic  distances.  Indeed  it  can  be  shown  that  the  next  per¬ 
turbation  is  of  order  e/c  . 

Substituting  (7.12)  into  (709)-  ,  we  obtain  t~,  along  the  crack 

line  x7  =  0  :  ^ 

“  00 
■ 

(7.13)  t23/t0  =  -  JjCOtl  -  C(S2  +  e2)"%]  cos(?x)  dc 

0 

From  (7.13),  it  can  be  deduced  that  as  e->“  (e  =  0)  t2,  gives  the 
classical  elasticity  solution.  For  e  =  0  (no  crack)  ^  t23  =  0  and 
for  e  ^  0  ,  there  is  no  singularity  along  the  x-axis.  Thez  non -dimen - 
sional  stress  T  =  t23/tQ  as  a  function  of  the  non-dimensional  distance 
X]/c  along  the  crack  line  is  shown  in  Fig.  4  for  two  different  values 
of  e(=  0. 1  and  0.01) .  From  these  figures,  it  is  clear  that  the  bound¬ 
ary  condition  at  the  crack  tip  is  satisfied  very  well,  especially  for 
e  =  0.01,  as  expected.  Outside  of  the  crack,  the  stress  reaches  a  maxi¬ 
mum,  afterwards  decaying  with  distance,  to  the  classical  elasticity 
solution.  In  our  previous  work  [3,5,11],  we  have  introduced  a  fracture 
criterium  based,  physically  realistic,  maximum  stress  hypothesis.  This 
gave  excellent  agreement  with  the  results  of  atomic  theories  and  experi¬ 
ment.  Further  study  of  this  problem  may  be  necessary  in  order  to  take 
account  of  the  inhomogeneity  caused  by  the  presence  of  the  crack. 


Figure  4>  stress  distribution 
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8  *  Cylindrical  Hole  (Defect) 


A  cylindrical  hole  of  radius  r  with  axis  x,  in  an  infinite 
elastic  medium  may  be  treated  as  a  two-dimensional  problem.  In  plane 
polar  coordinates, 


(8.1) 


=  r  cos  0  . 


With  the  Airy  stress  function  <}> (x^^) 

(8.2)  tn=  (J)^22  ,  X-22  =<^>H  » 


r2  =  r  sin  0 


*12  ""*,12 


equations  oi  equilibrium  are  satisfied  identically. 

The  strain  tensor  e^  must  obey  the  compatibility  equations 

ell,22  +  e22,ll  "  2  e12,12  =  0 


For  e^  ,  we  have  the  relations 
(8.4)  2 y  ek£  =  ~  v  , 


k,a  =  i»2 


and  for  a,,  from  (5 .11)  with  y  =  0  , 


(8.5) 


0-  -  e  V  )  tj^  -  ak£  , 


all  of  these  being  valid  in  rectangular  coordinates. 

Combining  (8.2)  to  (8.5),  we  obtain 
(8.6)  (1  -  e2V2)  V4  <j>  -  0 

But  from  (8.2)  to  (8.  5),  summing  over  k  and  l  ,  we  have 


(8.7) 

Consequently 

(8.8) 


Cl  -  eV)  V2*  *  V-u 


v  (V»u)  =  0 


whose  general  solution  is 


(8.9) 


7«u  =  C^  £nr  +  C2 


Since  V*u  can  be  replaced  by  its  expression  referred  to  polar  coordi¬ 
nates,  ~  ~ 

(8.10)  V*u  =  y  |p  (ru)  ,  u  =  u(r)er 


(8.9)  is  integrated  to  give 


(8.11) 


■u  =  -^r2(2  Jin  r  -  1)  +  j  C, 


r2  +  C 


19 


we  obtain 


Using  the  fact  that  u  must  vanish  as  r 

(8.12)  ru  =  C 
With  this,  V*u  =  0  and  (807)  are  integrated  to  give 

(8.13)  <j)  =  A-j^  +  A2  Jlnr  +  A3I0(r/e)  +  A^Or/e) 

where  A-,  to  A^  are  constants,  and  In  and  K  are  the  modified 
Bessel* s^funct ions.  To  calculate  the  stress  0  field  in  polar  coordi¬ 
nates,  we  observe  that 

(8.14)  trr  =  *n  =  7  af  »  t06  =  v  <j>  "  ^ 

The  stress  field  will  vanish  at  r  -  «>  if  A.  =  0  ,  and  it  will  have  the 
classical  elasticity  value  C  for  e  =  0  if5  A2  =  -2yC.  Hence 

trr  =  -  (2pC/e2p2)[l+ApK1(p)]  , 

(8.15)  t0e  =  (2yC/e2p2)[l4ApK1(p)  +  A  p2K0(p)], 


peu  =  C 

where 

(8.16)  p  =  r/e  ,  A  -  A4/2yC 

The  constant  C  is  connected  with  the  volume  change  at  r  =  0  ,  and 
constant  A  is  left  undetermined.  By  expanding  Kq  (p)  and  JC,  (p) 
near  p  =  0  ,  we  find  that  for-  A  =  -1  the  singularity  of  the  'stress 
field  can  be  reduced  from  p"z  to  £n  p  ,'t  p  =  0  ,  i.e.  in  the 
vicinity  of  p  =  0  ,  we  have 

t_  =  (uC/e2)  [)ln(p/2)  -  l  -  iKD)  , 

(8.17)  **  ^  2 

t00  *  (uC/e2)  f  to  (p/2)  +  KD] 

where 

-  #(1)  *  C  =  0.5772... 

C  being  the  Euler  constant.  In  terms  of  the  displacement  field,  these 
read 

trr  =  (ypu/e)  lfci(p/2)  -  j  -  4»(1)]  , 

(8.18) 

t00  =  (ypu/e)  (£n(p/2)  +  j  -  i|>(l)]  . 

On  physical  grounds ,  u  must  be  finite  at  r  =  0  ,  Since  interatomic 
distance  is  of  the  order  of  e  ,  we  see  tiiat  for  finite  u  ,  t  and 
t0fi  +0  as  p  -*-0. 
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'  wuewaj uuwi  «WW7CT_'ner«xfl3UWnim*jciUTsjri«  jrujnuru.  tuituj- 


Although  u  ■*■ »  as  p  ->-  0 ,  we  may  assume  that  the  lowest  value  of 
r  =  e  ,  i„e„  p  =  l  .  For  this  value  of  p  ,  both  t  _  and  taa  are 
finite.  rr  86 
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